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ABSTRACT 

We use the Chandra X-ray Observatory to study the dark matter halos of 34 mas- 
sive, dynamically relaxed galaxy clusters, spanning the redshift range 0.06 < z < 0.7. 
The observed dark matter and total mass (dark-plus-luminous matter) profiles can 
be approximated by the Navarro Frenk & White (hereafter NFW) model for cold 
dark matter (CDM) halos; for ~ 80 per cent of the clusters, the NFW model pro- 
vides a statistically acceptable fit. In contrast, the singular isothermal sphere model 
can, in almost every case, be completely ruled out. We observe a well-defined mass- 
concentration relation for the clusters with an intrinsic scatter in good agreement 
with the predictions from simulations. The slope of the mass-concentration relation, 
c oc M" ir /(1 + z) b with a — —0.45 ± 0.12 at 95 per cent confidence, is steeper than 
the value a ~ —0.1 predicted by CDM simulations for lower mass halos. With the 
slope a included as a free fit parameter, the redshift evolution of the concentration 
parameter, b = 0.71 ± 0.52 at 95 per cent confidence, is consistent with the same 
simulations (b ~ 1). Fixing a ~ —0.1 leads to an apparent evolution that is signifi- 
cantly slower, b = 0.30 ± 0.49, although the goodness of fit in this case is significantly 
worse. Using a generalized NFW model, we find the inner dark matter density slope, 
a, to be consistent with unity at 95 per cent confidence for the majority of clusters. 
Combining the results for all clusters for which the generalized NFW model provides 
a good description of the data, we measure a = 0.88 ± 0.29 at 95 per cent confidence, 
in agreement with CDM model predictions. 

Key words: cosmology: observations - dark matter - X-rays: galaxies: clusters 



1 INTRODUCTION 

One of the most remarkable results of cold dark mat- 
ter (CDM) simulations of structure formation is that 
the density profiles of dark matter halos on all resolv- 
able mass scales, from small satellites to the most mas- 
sive galaxy clusters, can be approximated by a universal 
profile, the so-called Navarro- Frenk- White (NFW) profile 
l|Navarro. Frenk fc Whitll 19951 . Il997h 
Po 



p(r) 



& 1 + t:) 2 



(1) 



Here r B is the characteristic scale radius of the halo and 
po is the central density. The model has a central, negative 
logarithmic density slope, referred to hereafter as the "inner 
slope" 



dlogp 



d logr 



= 1. 
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NFW also defined the concentration parameter, c, as the 
ratio of r2oo, the radius within which the matter density is 
200 times the critical density, and the scale radius: 

c=— . (3) 

' s 

They showed that smaller mass halos are more concentrated 
than the higher mass halos and interpreted this as reflecting 
the higher formation redshift of the lower mass systems. 

Copious numerical work has been devoted to 
testing these fundamental findings. To explain the 
mass-concentr a tion relation and its r edshi ft evolution, 
iBullock et all (|200ll ) and lEke et al.l i|200ll ) introduced 
simple, but highly success ful models for the f ormation of 
dark matter halos. In the IBullock et alj l|200l[ ) model, for 
example, only two parameters, K, which determines the 
initial concentration parameter of collapsing halos and F, 
the ratio of the initial collapse mass to the final virial mass 
of the halo at redshift zero, are required to approximately 
match the simulation predictions. More recently simple two- 
parameter power-laws have been used to characterize the 
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mass-concentration relation (|Dolag et al.ll2004l ; | Shaw et al.l 
l200rj ). 

When considering the inner slopes of dark matter halos, 
a useful generalization of the NFW formula is 



p(r) = 



Po 



(4) 



(e.g., lHernquisdll990l ; IZhao|[l996l ; |jing fc Sutoll2000l) . where 
the inner slope a is a free parameter. From their 



tial simulations, iMoore et al.l lrjl999T ) found a steeper cen- 
tral slope than NFW, with a = 1.5. For some years, the 
question of the precise value for the inner slope remained 
the topic of much debate. However, a consistent view has 
now emerged in which real dispersion i s expected between 
the inner slopes for i ndividual halos l|Klypin et alj l200ll ; 
iTasitsiomi et al.1 \2004 ) and where typical values for the 
inner slopes of clusters lie in the range a ~ 1.1 ± 0.4 
jNavarro et al.ll2004lD~iemand et ail2004l2005h . Recent an- 
alytical and numerical work to solve the Jeans equation sug- 
gests a value for the inner slope of a » 0.8 jAustin et al.l 
|2005| ; iDehnen fc McLaughlinll2005l ; lHansen fc Staddl2006j ). 

In summary, both a ti ght mass-concent ration relation 
(for recent results see e.g., IShaw et al.ll2006h and an inner 
density slope for dark matter halos in the range 0.7 < a < 
1.5 are central predictions of the CDM paradigm for struc- 
ture formation. 

In this paper we use Chandra X-ray Observatory data to 
study the mass profiles for 34 of the most massive, dynami- 
cally relaxed galaxy clusters known. Such clusters are among 
the most promising targets with which to check the central 
CDM predictions, being both dominated by dark matter 
(e.g., lAllen et al.l 12004 . and references therein) and hav- 
ing a size that allows us to resolve well within the scale ra- 
dius, even at high redshifts. We employ an analysis method 
which minimizes the need for priors associated with the use 
of parameterized models for the X-ray gas density and/or 
temperature profiles. The restriction to relaxed clusters also 
minimizes systematic effects associated with e.g. geometry 
and possible non-thermal pressure support. 

Previous X-ray studies have suggested broad agreement 
with the theoretically predicted mass-concentration rela- 
tion (|Pointecouteau et al.ll2005l ; IZhandl2006l ; IVikhlinin et al.1 



l200fj ; IVoigt fc FabianlbOOfj ), at least at low redshifts. Here, 
we increase the detail of this comparison and, for the first 
time, measure both the slope of the mass-concentration 
relation and its evolution. We show that tension may 
exist between the Chandra data and some simple mod- 
els based on CDM simulations for lower mass halos. 
With regard to the inner density slope, Chandra and 



good 


agreement with the CDM predictions (lAllen et al. 


2001 




Schmidt et al.1 l200ll; 


Allen et al.l |2002|; iLewis et al. 


2002 




20031; iBuote & Lewis 


2004; Arabadiis et al.l 2002; 



Arabadiis fc Bautj|2004l ; lAndersson fc Madeiskil 120041 ). Af- 
ter some initial controversy, there is also an emerg- 
ing consensus that strong gravitational lensing data sup- 
port inner dark matter density slopes of about unity 



in several cluster s (ISmith et al.l l200ll; IGavazzi et ail 



Sand et al.1 12002). |2004|; iBartelmann fc Meneghettil 
Itonl b 



2003 



2004 



Dalai fc Keetonl [20031 ; iMeneghetti et al.1 120051 ; ISand et al.l 
20051 ). We extend this work and show that the NFW model 



provides a good description of the total mass and dark mat- 



ter profiles for most relaxed clusters, rejecting completely 
the simple singular isothermal model (p oc 1/r 2 ). We also 
obtain a robust result on the inner slope for the ensemble of 
clusters. 

A flat ACDM reference cosmology is assumed with 
Hubble constant Ho = 70kms _1 Mpc _1 and matter density 
fl m = 0.3. In a few cases cluster masses are quoted with a 
Hubble constant scale h — i/o/100kms _1 Mpc -1 . 



2 METHOD 

2.1 Target selection 

Our target clusters are the most massive, dynamically re- 
laxed clusters known in the redshift range < z < 0.7. 
They form a restricted set of the clusters used by Allen 
et al. (2007, in preparation) to study the evolution of the 
X-ray gas mass fraction and constrain cosmological param- 
eters. In detail, we have used only the targets from that 
study for which the temperature is measured in at least four 
independent bins, which permits reliable measurements on 
the inner density slopes. The target clusters all have mass 
weighted X-ray temperatures measured within the radius 
7"25<xEl kTzsoo 5 keV (Allen et al. 2007, in preparation). 

The clusters exhibit a high degree of dynamical relax- 
ation in their Chandra images, with sharp central X-ray 
surface brightness peaks, regular elliptical X-ray isophotes 
and minimal isophote centroid variations. The clusters show 
minimal evidence for departures from hydrostatic equilib- 
rium in X-ray pressure maps (Million et al. , in prepa- 
ration). The exceptions to this are RXJ1347. 5-1145, and 
MACSJ0744.9+3927, for which clear substructure is ob- 
served between position angles of 90-190 degrees and 210- 
330 degrees, respectively. These regions, associated with ob- 
vious substructure, have been excluded from the analysis. 
The restriction to clusters with the highest possible degree 
of dynamical relaxation (for which the assumption of hydro- 
static equilibrium should be most valid) minimizes system- 
atic scatter and allows the most precise test of the CDM 
model predictions (e.g., iNagai et al.ll2007h . 



2.2 Observations, data reduction and spectral 
analysis 

The Chandra observations were carried out using the Ad- 
vanced CCD Imaging Spectrometer (ACIS) between 1999 
August 30 and 2005 June 28. The standard level-1 event 
lists produced by the Chandra pipeline processing were re- 
processed using the CI AO (version 3.2.2) software package, 
including the latest gain maps and calibration products. Bad 
pixels were removed and standard grade selections applied. 
Where possible, the extra information available in VFAINT 
mode was used to improve the rejection of cosmic ray events. 
The data were cleaned to remove periods of anomalously 
high background using the standard energy ranges and time 
bins recommended by the Chandra X-ray Center. The net 
exposure times after cleaning are summarized in Table [1] 



1 ^2500 is the radius within which the mean mass density is 2500 
times the critical density of the Universe at the redshift of the 
cluster. 
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Table 1. Summary of the Chandra observations. Columns list the target name, redshift, observation date, detector used, observation 
mode, net exposure after all cleaning and screening processes, and coordinates from the X-ray centres used in the analysis. Where multiple 
observations of a single target have been used, these are listed separately. Redshifts for the MACS clusters are from Ebeling et at. 2006, 
in preparation, and will appear in full in the published article. 
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The d ata have b e en an alysed using techniques dis- 
cussed by lAllen et all (|2004T ) and references therein. In 
brief, concentric annular spectra were extracted from 
the cleaned event lists, centred on the coordinates 
listed in Table [1] Emission associated with X-ray point 
sources or obvious substructure (Section 12. ip was ex- 
cluded. T he spect r a we re analysed using XSPEC (ver- 
sion 11.3: lArnaudl Il996h, t he MEKAL plasma emission 
code dKaastra fc Mewej 19931 ; incorporating the Fe-L calcu- 



lations of lLiedahl et alJliMMai and tne pnotoeiectric absorp- 
tion models of Balucinska-Church fc McCammonl l|l992h . 
We have included standard correction factors to account for 
time- dependent contamination along the instrument light 
path. In addition, we have incorporated a small correction to 
the High Resolution Mirror Assembly model in CIAO 3.2.2, 
which takes the form of an 'inverse' edge with an energy, 
E=2.08keV and optical depth r = -0.1 (H. Marshall, pri- 
vate communication) and boosted the overall effective area 
by six per cent, to better match later calibration data (A. 
Vikhlinin, private communication). Only data in the 0.8 — 7.0 
keV energy range were used in the analysis (the exceptions 
being the earliest observations of 3C 295, Abell 1835 and 
Abell 2029 where a wider 0.6 to 7.0 keV band was used). 



For the nearer clusters (z < 0.3), background spectra 
were extracted from the blank-field data sets available from 
the Chandra X-ray Center. These were cleaned in an iden- 
tical manner to the target observations. In each case, the 
normalizations of the background files were scaled to match 
the count rates in the target observations measured in the 
9.5-12keV band. Where required, e.g. due to the presence of 
strong excess soft emission in the field of Abell 2029, a spec- 
tral model for any unusual soft background emission was in- 
cluded in the analysis. For the more distant systems (as well 
as for the first observation of Abell 1835, the ACIS-I obser- 
vation of Abell 383, and the observations of Abell 2537, RXJ 
2129.6+0005 and Zwicky 3146) background spectra were ex- 
tracted from appropriate, source free regions of the target 
data sets. (We have confirmed that similar results are ob- 
tained using the blank-field background data sets.) In order 
to minimize systematic uncertainties and due to the specific 
goals of this work, we have restricted our spectral analysis 
to radii within which systematic uncertainties in the back- 
ground subtraction (established by the comparison of differ- 
ent background subtraction methods) are smaller than the 
statistical uncertainties in the results. All results are drawn 
from spectral analyses limited to ACIS chips 0,1,2,3 and 7 
which have the most accurate calibration, although ACIS 
chip 5 was also used to study the soft X-ray background 
in ACIS-S observations. We do not attempt to extend our 
analyses to larger radii using the data from other chips. 



2.3 Cluster mass profile measurements 

2. 3. 1 Basics of the mass analysis 

Under the assumptions of hydrostatic equilibrium and 
spherical symmetry, the observed X-ray surface brightness 
profile and the deprojected X-ray gas temperature pro- 
file may together be used to determine the X-ray emit- 
ting gas mass and total mass profile of a galaxy cluster. 
For this analysis, we have used an enhanced version of 
the Cambridge X-ray d e proje ction code described by e.g. 
IWhite. Jones fc Formanll|l997h . This method is particularly 
well suited to the present study in that it does not require 
approximate fitting functions for the X-ray temperature, gas 
density or surface brightness when measuring the total, grav- 
itating masfl The use of such functions introduces priors 
into an analysis which can complicate the interpretation of 
results and, in particular, the estimation of measurement 
errors. 

We have carried out two separate mass analyses: firstly 
(method 1) an analysis in which the total mass profile (dark 
plus luminous matter) was modelled using either an NFW 
or singular isothermal sphere. Detailed results on the X- 
ray emitting gas profiles were also determined at this stage. 
Secondly (method 2) an analysis in which the total mass 
profile was modelled as the sum of three parts: the dark 
matter halo (fitted with a generalized NFW profile), the X- 
ray emitting gaseous halo (approximated, for the purposes 
of this analysis only, with a beta model), and the optically 
luminous mass of the cD galaxy (approximated with a Jaffe 
or de Vaucouleurs model). 

The normalization, po, of the generalized NFW mass 
model (eq. [3} is usually written as 



p0 = Pcrit 8 C 

where 

p clit = 3H{z) 2 /8nG 



(•») 



(6) 



is the critical density at the redshift z of the cluster. G is 
the gravitational constant and the Hubble parameter H(z) 
in the flat ACDM reference cosmology is defined by 



H{zf = H 2 ((1 + z) 3 fi m + 1 - fi n 



(7) 



where f2 m is the matter density in units of the critical den- 
sity. 

The amplitude S c depends only on the concentration pa- 
rameter c = -^ a - (r2oo is the radius within which the average 

*"20g v to 
mass density is 200 times p C rit)- For convenient computation, 

8 C can be written using the Gauss hype rgeometric function 
F(a,b;c;z) jAbramowitz fc Stegu5ll965h 



3 — a 



F(3-a,3-a;4-a;-y) 



(8) 



Separate photon-weighted response matrices and effec- 
tive area files were constructed for each region using cali- 
bration files appropriate for the period of observations. The 
spectra for all annuli for a given cluster were modelled simul- 
taneously in order to determine the deprojected X-ray gas 
temperature and metallicity profiles, under the assumption 
of spherical symmetry. 



5c 



200 c 3 



(9) 



(|Wvithe et alfcOOll ). For the original NFW mass model with 

2 As discussed in the text, priors are required when modelling the 
dark matter and X-ray and optically luminous matter components 
separately 
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q = 1, one obtains <fr(c) = ln(l+c) — c/(l+c) (jNavarro et al.l 
1 19961 ). 

With mass analysis method 2, the component of 
the total mass distribution representing the X-ray emit- 
ting gas component was desc ribed by a beta-model 
|Cavaliere fc Fusco-Femianoiri978l ) 



,{r) = po.g 



1 + 



(10) 



where po.gas is the central gas density, r Cigas is the core ra- 
dius of the gas profile and /3 is the slope parameter. (We 
stress that the detailed results on the X-ray emitting gas 
profiles, used for example in the measurement of cluster gas 
mass fractions, are determined with method 1 and involve 
no assumption about the functional form of the gas profile. 
The use of the beta-model with method 2 simply approxi- 
mates the contribution of this mass component to the total 
mass, allowing us to recover the dark matter profiles. Al- 
though the beta model does not provide a precise match to 
the X-ray gas mass distribution in all cases, the systematic 
uncertainties in the dark matter profiles that result from its 
use are small, primarily because the X-ray emitting gas con- 
tributes only ~ 12 per cent of the total mass; Allen et al. 
2007, in preparation). 

All of the target clusters have a single, optically domi- 
nant galaxy at their centres. With method 2, we accounted 
for the mass of stars in the cent ral galaxy using a lJaffd (|l983l ) 
model i|Sand et al.ll2002l . 120041 ). This was added to the gen- 
eralized NFW potential for the dark matter halo and the 
beta-model for t he X-ray em itting gas to obtain the total 
mass profile. The IJaffel (|l983l ) model is 



PA 



P0,3 



(7T) 2 (l + 77) 



(11) 



with central density po,J and core radius r c . The total mass 
for the model is finite, Mj = 4nr^ po,j- Since including this 
model into the analysis has only a small effect on the re- 
sults (Section I4.2|l but requires two parameters, we chose 
to fix these two parameters to sensible values. We adopt 
R e = 0.76 r c = 25 kp c as a typical effective radius (e.g., 
ISand et alj I2002L 120051) . For th e well-studied cD galaxy in 
MS2137-2353. ISand et ail |2002l ) measure a total V-band lu- 
m inosity of Ly = 4.16 x 10 11 L©. Using the M/Ly relation 
of iFukugita et al l (|l998l ) 

M 

— = 4.0 + 0.38 (i g - lOGyr), (12) 
Lv 

where i g is the age of the galaxy (we assume a formation red- 
shift Si — 2.0), we estimate the total stellar mass associated 
with the cD galaxy of 1.14 x 10 12 M Q . The central galaxy 
for each cluster in the sample is assumed to have this stellar 
mass. 

We hav e also carried out a repeat analysis in which the 
IJaffel i| 19831 ) model for the stellar mass associated with the 
centra l, dominant galaxies was replaced by a ide Vaucouleursl 
( 1948) model. This analytical model is described by the sur- 
face mass density 



= E c e 



-7.67[(r/flo) 1/4 -l] 



(13) 



where E c is the surface density at the effective (or half- 
mass) radius R e . The 3-dimensional mass profile of the de 



Table 2. Regions with residual substructure that were down- 
weighted in the mass analysis. A systematic uncertainty of ±30 
per cent has been added in quadrature to all temperature mea- 
surements made within radii -R su b- 

Rsub (kpc) 



Abell 1795 


73.3 


Abell 2029 


30.0 


Abell 478 


14.6 


PKS0745-191 


53.0 


Abell 1413 


38 


Abell 2204 


76.8 


Abell 383 


38.7 


RXJ1504. 1-0248 


79.1 


RXJ2129.6+0005 


40.4 


Zwicky 3146 


242 


Abell 2537 


41.1 


MACSJ2229. 8-2756 


42.0 


MACSJ0947.2+7623 


40.1 


MACSJ1931.8-2635 


41.5 


MACSJ1115. 8+0129 


83.4 


MACSJ1532.9+3021 


42.3 


MACSJ1621.6+3810 


43.1 



Vaucouleurs model was calculated bv lYounsl l)l976J ) and has 
a shallower central slope close to unity. However, similar 
results on the cluster dark matter profiles were obtained in 
all cases, showing that the precise choice of the galaxy model 
has a negligible effect on the results. 



2.3.2 Best-fitting values and confidence limits 

Given the observed surface brightness profile and a particu- 
lar parameterized model for the total mass, the deprojection 
code is used to predict the temperature profile of the X-ray 
gas. This model temperature profile is compared with the 
observed spectral, deprojected temperature profile and the 
goodness of fit is calculated using the sum over all temper- 
ature bins 



E 

all bins 



Tobs Tniod 
v CT obs 



(14) 



where T bs is the observed, spectral deprojected tempera- 
ture profile and T mo dd is the model, rebinned to the same 
spatial scale using flux weighting^ 

For each mass model, we determine the best-fitting pa- 
rameter values and uncertainties via x 2 minim ization. We 
use t he LEASQR Levenberg-Marquardt routine (Marquardt 
I1963T ) from the GNU Octave Repository, available online at 
http://sourceforge.net. (A full Monte Carlo analysis, as de- 
scribed above, is run for each set of parameter values exam- 
ined in the minimization procedure.) 



3 In detail, we use the median model temperature profile deter- 
mined from 100 Monte-Carlo simulations. The outermost pres- 
sure, at the limit of the X-ray surface brightness profile, is fixed 
using an iterative method that ensures a smooth, power law pres- 
sure gradient in these regions. The model temperature profiles, for 
radii spanned by the spectral data, are not sensitive to reasonable 
choices for the outer pressures. 
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For a number of the clusters, the Chandra images 
indicate small levels of residual substructure in the in- 
nermost regions, which probably result from 'sloshing' 
of the X- ray emitting gas within the central poten- 



tials ( e.g., Markevitch et af]|200ll ; lAscasibar fc Markevitchl 



[20061 : lAllen et alJ Il992h and/or interactions between cen- 
tral radio sources and the surrounding intraclust e r gas 



Birza n et al 
20051 : 



2004 



Dunn et al 



(e.g., Bohring er et al.lll9 93: Fabi an et all 12 000. 2003, 200 



Dunn & Fabian 2004: Forman et al 



20051 : iRaffertv et alJ 120061 ; Allen et al 



2006). The regions affected by such substructure are listed 
in Table [5] A systematic uncertainty of ±30 per cent has 
been added in quadrature to the spectral results from these 
regions which leads to them having little weight in the mass 
analysis. 



3 RESULTS 

3.1 Total mass profiles: NFW versus singular 
isothermal sphere 

Table [3] summarizes the results from the initial mass analy- 
sis (method 1) in which the total mass profiles (dark-plus- 
luminous matter) were modelled using either an NFW (in- 
ner slope a = 1 fixed) or singular isothermal sphere model 
(p(r) = A/r 2 , with the normalization A free). We see that 
for most clusters the NFW model provides a reasonable de- 
scription of the total mass profiles. For 27 out of 34 clusters, 
the probability of the y 2 / DOF value based on a \ 2 dis tribu- 
tion is 0.05 or better (e.g.. lBevington fc Robinsonlll992h and 
only for Abell 383 and 1413 does the probability drop below 
0.001. In contrast, the singular isothermal sphere model can 
be firmly rejected for most clusters in the sample. 

Combining the results for all clusters, the NFW model 
gives a total x 2 of 171.0 for 95 de grees of freedom (DOF) 
(with Abell 383 and 1413 contributing a total x 2 of 47.9). 
The singular isothermal sphere gives \ 2 — 3031.1 for 129 
DOF, indicating an extremely low model probability. 



3.2 Dark matter profiles: The NFW model and 
mass-concentration relation 

Table [4] summarizes the results from the analysis with 
method 2, in which the cluster mass distributions were sep- 
arated into three parts: the dark matter halo (fitted with a 
generalized NFW profile), the X-ray emitting gaseous halo 
(approximated with a beta model), and the optically lumi- 
nous mass of the central dominant galaxy (approximated 
with a Jaffe or a de Vaucouleurs model). In the first case, 
we examined models in which the inner slope of the dark 
matter profile was fixed at a = 1 i.e. the standard NFW 
model. 

Our first conclusion is that, as with the analysis of the 
total mass profiles, the NFW model provides a good overall 
description of the dark matter profiles in the clusters. For 
27 out of 34 clusters, the x 2 /DOF value has a probability 
of 0.05 or better and only for Abell 383 and 1413 does the 
model fail significantly. Combining the results for all clus- 
ters, the NFW model gives a total x 2 of 176.1 for 95 de grccs 
of freedom (DOF) (with Abell 383 and 1413 contributing a 
total x 2 of 48.0; note that the inclusion of the separate mass 




10 15 
M vir (h-iM ) 



Figure 1. The virial mass-concentration relation for the dark 
matter halos. The solid line shows the best-fitting model with 
c = 7.55 ± 0.90, a = -0.45 ± 0.12 and b = 0.71 ± 0.52 (95 
per cent confidence limits). T he dashed line show s the prediction 
from the CDM simulations of lShaw etafl <200d) . with c = 6.47, 
a = —0.12 and b = 1.0 fixed. The two clusters for which the NFW 
model fails to provide a reasonable description of the Chandra 
data, Abell 383 and 1413, are plotted with open square symbols 
and have been excluded from the fits. 



components for the X-ray emitting gas and stars introduces 
no additional free parameters in the fits.) The concentration 
parameters, c for the dark matter profiles are slightly lower, 
and the scale radii slightly higher, than for the total mat- 
ter distributions. However, in general the results are quite 
similar. 

As discussed in Section [T] one of the central predictions 
from simulations of CDM halos is the mass-concentration 
relation. In order to allow for the most direct comparison 
with theory, we have used the values in Tables [3] and [4] to 
calculate the concentration parameters c v ir and virial masses 
Mvir measured within the virial radii for the clusters. We 
adopt the d e finitio ns of virial radius and virial mass used by 
IShaw et al l (|2006T ) 



4 3 

M v ir = -TTrVir A c (z) pcrit(z), 



(15) 



where the virial ov erdensity is given by A c = 178f2 m (z) ' 45 
([Lahav et al.lll99l[ ). 

Both the virial radii and virial masses are calculated for 
the total mass model, including all mass components. The 
concentration parameters are defined by the ratios of these 
virial radii and the scale radii of the dark matter compo- 

Table [5] lists the results on Cvir 



nents, c v i r 



total / dark 



and Af v ir for the clusters using the standard NFW model 
(« = !)• 

Fig. [T] shows the variation of c v i r with M v i r measured 
from the Chandra data. A tight relation is observed with 
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Table 3. Results from fits to the total mass profiles with NFW (slope a = 1 fixed) and singular-isothermal (SI) models (analysis method 
1). Columns list the number of temperature bins used in the analysis of each cluster, the scale radius r s (in kpc) and concentration 
parameter c = r20o/f*s for the NFW models, and the goodness of fit (x 2 /DOF) for both models. The SI model provides a poor description 
of the data. Quoted uncertainties are 68 per cent (Ax 2 = 1.0) confidence limits. 



TOTAL MASS (LUMINOUS PLUS DARK MATTER) 
Number of NFW MODEL SI MODEL 

kT bins r s c X 2 /DOF X 2 /DOF 



Abell 1795 


5 




4 5Q+ - 86 
-0.62 


2.06/3 


49.2/4 


Abell 2029 


6 


0.28±° f 2 


6 Q5+ - 30 
°- a,J -0.39 


4.73/4 


303/5 


Abell 478 


7 


U - OU -0.06 


4 48+ 0,40 
^•*°-0.25 


8.85/5 


615/6 


PKS0745-191 


7 






5.13/5 


42.1/6 


Abell 1413 


6 




A on+0.69 


24.6/5 


91.5/6 


Abell 2204 


5 




9.84+j;? 4 . 


4.38/3 


16.6/4 


Abell 383 


5 


0.48+0;} 2 


q 7^+0-71 
a - 10 -0.48 


23.3/4 


90.0/5 


Abell 963 


5 




A 73+ -8 4 


6.16/3 


97.8/4 


RXJ0439.0+0521 


4 


u ' ly -0.04 


'• DM -1.24 


0.30/2 


28.3/3 


RXJ 1504. 1-0248 


5 


0.54lJ*j 


4 3S+ 1 - 06 


0.91/3 


28.1/4 


RXJ2129.6+0005 


5 


0.38t o f 7 


4 c iQ+ 2 ' 18 
*' oa -1.84 


1.53/3 


14.7/4 


Abell 1835 


5 




4 20 +0 ' 43 
^• zu -0.30 


12.5/3 


302/4 


Abell 611 


5 


_i_n on 
U " r>Z -0.10 


o_i fin 
5.39j:J;» 


2.35/3 


25.7/4 


Zwicky 3146 


5 


0.99-0.57 


z -' i -2.71 


0.79/3 


11.4/4 


Abell 2537 


5 




4-86t 2 '° 6 2 


3.42/3 


22.2/4 


MS2137.3-2353 


6 


o.islH 2 




5.52/4 


235/5 


MACSJ0242. 6-2132 


1 


. , in nr; 

1 o+ u - UD 
u - la -0.03 


' - oy -1.04 


2.92/2 


62.1/3 


MACSJ1427.6-2521 


4 


n i c+0.08 
U - 1O -0.05 


Q 1 q + 2.28 
°- 18 -1.92 


4.77/2 


18.1/3 


MACSJ2229. 8-2756 


4 


U - 1O -0.05 


8.43+J; 27 


1.69/2 


6.8/3 


MACSJ0947.2+7623 


1 


n 32+ - 18 

U - JZ -0.11 


e-oill:!! 


2.53/2 


23.6/3 


M ACS J1931. 8-2635 


4 


5i+ - 95 

u ' d -0.20 


4 05 +i.54 

^■ ud -1.91 


1.98/2 


33.6/3 


MACSJ1115. 8+0129 


4 


1.35_o.9o 


2.l5l 2 ;« 


1.23/2 


13.7/3 


MACSJ1532. 9+3021 


5 


34+ 16 
u -° -0.09 


5.29li' 22 


1.76/3 


41.8/4 


MACSJ0011. 7-1523 


4 


49+ 17 


+0.79 
°- '-0.63 


6.75/2 


168/3 


MACSJ1720. 3+3536 


4 


u -°°-0.09 


k ni+1-06 
°- zl -0.98 


1.19/2 


49.8/3 


MACSJ0429. 6-0253 


4 


16+ 004 
u ' 1D -0.03 


8.50+i; 4 ° 


1.84/2 


25.8/3 


MACSJ0159. 8-0849 


4 


3 7 +0.i0 
u -°'-o.io 


r- 9 o + 1.19 
tl -^-0.74 


4.41/2 


52.9/3 


MACSJ0329. 7-0212 


4 


29+ - 07 
u,z3 -0.05 


c z.o+0.76 
5-48_ 67 


6.24/2 


86.0/3 


RXJ1347.5-1144 


6 


45+ 010 
u '^°-0.07 


5 fi3+ - 54 


16.1/4 


257/5 


3C295 


4 


UiiO -0.02 


°- od _0.89 


1.12/2 


61.3/3 


MACSJ1621. 6+3810 


1 


25+ ' 17 
u ' zo -0.10 


°-°'-1.91 


4.00/2 


17.6/3 


MACSJ1311.0-0311 


4 


31+ 015 
u -° -0.09 


4 91+1-21 
*- al -1.10 


1.54/2 


25.4/3 


MACSJ1423.8+2404 


4 


17 +0.03 
"•^'-0.02 


o 97 +0.72 


0.05/2 


74.1/3 


MACSJ0744.9+3927 


4 


32+0-18 
u -° -0.09 


4.88li;« 


4.40/2 


40.9/3 



a clear trend for decreasing concentration parameter with 
increasing mass. 

The literature contains a variety of simple analyt- 
ical forms to describe the expected form of the mass- 
concentration relation, based on CDM numerical sim- 
ulations. The sim plest of these is a power-law model 
l|Dolag et al.ll2004l ) 

r ■s c f M vir V 
Cvir(z) = m j . (16) 

Sh aw et all (|2006t ) present one of the largest statistical stud- 
ies of cluster-sized CDM halos to date. Fitting a z = 0.05 
snapshot of their simulated data with the lDolag et al.l (|2004l ) 
power law model, these authors find Co = 6.47 ± 0.03 and 
a — —0.12 + 0.03 (68 per cent confidence limits). The dashed 
line in Fig. [1] shows the model defined by these parameters, 
overlaid on the Chandra results. The model clearly provides 



a poor fit to the observations, \ 2 = 125.9 for 30 DOF, lying 
systematically below the data at lower masses and above it 
in the highest mass range. Note that the two clusters for 
which the NFW model provides a formally unacceptable de- 
scription of the Chandra data, Abell 383 and 1413, have 
been excluded from the fit and are plotted with open square 
symbols. N o te als o that the mass-concentration relation of 
[ Shaw et al l (120061) is in good agreement with t he models of 
[ Bullock et al.l d200ll) (with the parameters from lDolag et all 
l2004h and lEke et al.1 |200ll ). 

One should remember that the Chandra observations 
reported here are for the most massive, dynamically re- 
laxed clusters known within z < 0.7. (Our targets repre- 
sent the most relaxed ~ 20 per cent of clusters in this 
mass and redshift range.) We have, therefore, also deter- 
mined the best-fitting power-law parameters appropriate 
for the most relaxed 20 per cent of simulated halos with 
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Table 4. Results on the dark matter profiles for the clusters, from the fits including separate components for the X-ray emitting gas and 
dominant cluster galaxy (analysis method 2). For the standard NFW model (inner slope a = 1 fixed) we list the scale radius r B (in kpc) 
and concentration parameter c = r20o/r B - For the generalized NFW model, we give results on the inner slope a. For those cases where 
the formal best-fitting values for the slope are negative, we only list the 95 per cent upper limit. Where confidence limits are absent 
in the table, the fit results were consistent with a = 0. For both models, we list the results on the goodness of fit (x 2 /DOF). Quoted 
uncertainties are 68 per cent (Ax 2 = f-0) confidence limits. 



DARK MATTER PROFILES 
NFW (a = 1) GENERALIZED NFW 





r s 


C 


X 2 /DOF 


a 


x 2 /DOF 


A Kail ~\7QX> 
.MUCH 1 ( yj 


41+ 013 
U - 41 -0.09 


A 4 r- + 0.85 
4 - 40 _0 .77 


2 21/3 


s 1 

< l.UO 


1 ^7 /9 
l.JI J z 


Abell 2029 


28+ ' 03 


6 63+ ' 34 

o - uo -0.37 


4.47/4 


1 16+ 022 

1 ' lo -0.34 


4.31/3 


Aholl 478 
r\. [Jell *± 1 O 


c ifi+ 09 
u-oo_ og 


q qi+0 36 
- a -0.33 


QO 

u.yu/ 


1 10+ 012 

i - iU -0.21 


fi fifi /4 


PKSf-^'S 1Q1 


3fi+ 013 


5 85+ 1 ' 55 


5.26/5 


0.78+0.64 


5.17 /4 


A Kail 1413 
_-_DCll ±_;±0 


n 43+ 014 

u '*J-0.09 


4 4 o+0.78 
*-*°-0.75 


Z^_. 0/ 




99 n /zl 

ZZ.U/ 4_ 


A Kail 9904 
iADCll ZZU^i 
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Mvir > 2 x 10 14 /i _1 M Q ) in the IShaw ei~ai1 l|200fj ) study. 
(In detail, this is done by selecting only those clusters 
with a substructure fraction, as defined by those authors, 
f s < 0.045). The resulting power- law fit parameters are 
Co = 6.8 and a = —0.16. Although giving slightly better 
agreement with the observations, this model still provides a 
poor description of the Chandra data with x 2 = 102.9 for 
30 DOF. 

Physical motivation for the power-law model and its 
(l_fz) _1 redshift dependence is presented bv lBullock et all 
(|200ll ). In their model they assume that the ratio of the 
virial mass, M v j r , to the volume inside the scale radius, V s , is 
proportional to the matter density p m (2coii) at the redshift, 
Zcoii, when the cluster formed: 



. . Mvir Pm T vil 3 , . 

Pm(-Zcoll) OC —rj— OC 3 = p m Cvjr. (17) 

Since the matter density p m was larger in the past, p m oc 
(1 + z) 3 , the typical concentration parameter for a halo of 
mass Myjr is s malle r at higher redshifts, c v i r oc (1 + z) . 
iBullock et alj (|200lf ) associate each redshift with a typical 
collapsing mass (defined as a fraction F of the final halo 
mass) through a critical value of the variance a(z,M) of 
the density fluctuations. Since a(M) for a ACDM model 
is a power-law for masses M v i r < 10 13 Mq/F, this model 
yields a power-law for the mass-concentratio n relation (up 
to thi s mass scale) with two free parameters. IBullock et alj 
l|200ll ) suggest F = 0.01 (although F = 0.001 is also 
used in the more recent literature). At the high-mass end, 
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Figure 2. The measured concentration parameters for the dark 
matter halos as a function of redshift. The solid line shows the 
best-fitting model with c = 7.55 ±0.90, a = -0.45 ±0.12 and b = 
0.71±0.52 (95 per cent confidence limits). Th e dashed line shows 
the prediction from the CDM simulations of Shaw ct al. d2006h 
with Co = 6.47, a = —0.12 and b = 1.0. The apparent absence 
of redshift evolution in the data may be an artifact induced by 
the steep slope of the underlying mass-concentration relation (see 
text for details). Symbols as in Fig. [T] 



Mvir > 10 13 Mq/F, the model predicts a change in the 
power-law slope but, unfortunately, also becomes unrealistic 
as the linear evolution of density fluctuations stalls in the 
ACDM model, preventing high-mass clu sters from forming 
jBullock et al.ll200ll; iKuhlen et al.ll2005l ). 

Although the Bullock et al.l " i|200lf ) model and the power 
law model of iDolag et al.l (|2004j) are attractive in terms of 
their simplicity, the limitations of such models for describing 
the detailed properties of cluster halos should not be over- 
looked. Firstly, real clusters contain X-ray emitting gas and 
stars, as well as dark matter. CDM-only simulations do not 
include cooling and feedback processes which affect the bary- 
onic mass components and modify the overall mass distribu- 
tions. Secondly, to date, studies of the mass-concentration 
relation have included very few halos a t the larges t mass 
range spanned by the Chandra data; e.g. JShaw et al.l (|2006l ) 
have only a single cluster with M v i r > 10 15 fo -1 Mq. Finally, 
Zhao et al] (|2003l) a rgue that at the highest masses, the 
Bullock et al.1 (|200lT ) model may over-predict the evolution 
of c v i r with redshift. 

Motivated by such considerations, we hav e introduced 
additi onal freedom into the power law model of IDolag et al.l 
l|2004l ). As well as having co and a as free fit parameters, we 
also allow the redshift evolution to evolve as (1 + z)~ b , with 
b free. 



CO 



Mvir 



(1 + z) b V 8 x 10 14 h~ 1 M Q 



(18) 



The results from a fit with this model, with Co, a and b all 



free, are shown by the solid line in Fig. [T] The model pro- 
vides a significantly improved description of the data with 
X 2 = 41.5 for 29 degrees of freedom and best-fitting param- 
eters co = 7.55 ± 0.90, a = -0.45 ± 0.12 and b = 0.71 ± 0.52 
(95 per cent confidence limits; similar results are obtained 
from a Monte Carlo analysis of the (c v i r |M v i r ) data wherein 
the results for individual clusters are scattered according to 
their measurement errors). Although the normalization Co 
at 8 x 10 14 fo ^Mp) and the red shift evolution, b, are consis- 
tent w ith the lShaw et al.l (|2006l ) simulations and lDolag et al.l 
l|2004h model, the mass-concentration relation is noticeably 
steeper. 

Fig. [2] shows the variation of concentratio n parameter 
with redshift. At low redshifts, the fit to th e Shaw et al.l 
l|2006l) simulated data with the IDolag et all (|2004l ) model 
(co = 6.47, a — —0.12, 6=1) provides a reasonable 
ma tch to the data (in agreement with the conclusions drawn 
by iPointecouteau et all 120051 ; IZhand 120061 ; IVikhlinin etafl 
l2006t IVoiet fc Fabianll2006l ). At higher redsh i fts, h owever, 
the observed c vlI values exceed the lShaw et al.l (|2006h model 
predictions. The most striking feature of Fig. [5] is an appar- 
ent absence of redshift evolution in the concentration pa- 
rameter. In d eed, w hen fixing the slope parameter a to the 
Sh aw et al] (|2006l ) value of a — —0.12, a fit with Co and 
b free gives a redshift evolution parameter b = 0.31 ± 0.49 
(95 per cent confidence limits), consistent with no evolution 
or even positive evolution of c v i r with M V [ T . However, note 
that the x 2 f° r this fit is significantly worse (X 2 = 74.7 for 
30 DOF; Ax 2 = 33.2) than for the fit in which the slope, a 
of the mass-concentration relation is also included as a free 
parameter (see above). We conclude that any analysis of 
redshift evolution in the mass-concentration relation should 
(at least) explore the full parameter space discussed above, 
with Co, a and b included as free parameters. 

Finally, we have estimated the systematic scatter that 
may be present in the observed mass-concentration relation. 
This was carried out by modifying the \ 2 estimator to in- 
clude an additional systematic uncertainty, and increasing 
the size of this systematic uncertainty until the reduced \ 2 
value became equal to unity. Based on the fit with Co , a and b 
included as free parameters, we estimate an intrinsic system- 
atic scatter in the data of Aiog(c) ~ 0.1, in g ood agreement 
with the predictions from simulations (e.g., iBullock et al.l 
|2001| ; IWechsler et~ai1l2002l ) . 

In summary, the key results on the mass-concentration 
relation are 1) the presence of a tight, observed mass- 
concentration relation for massive, dynamically relaxed clus- 
ters. The estimated intrinsic, systematic scatter in this re- 
lation is consistent with the predictions from CDM simula- 
tions. 2) We observe a slope a = —0.45 ± 0.12 (at 95 per 
cent confidence) for the mass-concentration relation that is 
significantly steeper than pre dicted by CDM si mulations for 
lower mass halos (a ~ -0.12; IShaw et al.ll2006] ). 3) The red- 
shift evolution of the observed mass-concentration relation, 
b = 0.71 ± 0.52 at 95 per cent confiden c e, is c onsistent with 
the value of b = 1 in the lBullock et al.1 l|200ll ) model. 



3.3 The inner slope of dark matter profiles: 
generalized NFW models 

For the final stage of our analysis, we have included the in- 
ner slope of the dark matter profile (a, in the generalized 
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Table 5. Results on the concentration parameters c v ; r = 
r vi r tal /^ ark and viria l masses M vir (in units of 10 14 h' 1 M Q ) for 
the standard NFW model (inner slope a = 1 fixed). Uncertainties 
are 68 per cent confidence limits. If no error bar is quoted, the 
upper limit is unbounded and the lower limit is zero. 
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NFW model) as an additional free parameter in the fits with 
method 2. (Separate components to model the mass contri- 
butions from the X-ray emitting gas and the optically lumi- 
nous matter in the central, dominant galaxies were included 
in the fits.) 

The results on the inner mass profiles are summarized in 
Table|4] In terms of the goodness of fit, 27 out of 34 clusters 
give x 2 /DOF values with a probability of 0.05 or better. 
Of these 27 clusters, 21 have an inner slope consistent with 
unity at 68 per cent confidence. 

In order to obtain a combined result on the inner dark 
matter profile, we have summed together the results on \ 2 
as a function of a for the 27 clusters for which the gener- 
alized NFW model provides a reasonable description of the 
data. The results are shown in Fig. [3] (Note that we have 
subtracted the overall minimum summed x 2 , of 75.37 for 72 
DOF; 127 temperature measurements, 2 free parameters for 
27 clusters and one free slope parameter). 
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0.5 



Figure 3. The summed \ 2 values as a function of the inner den- 
sity slope a for the 27 clusters for which the generalized NFW 
model provides a reasonable fit to the data. The dashed line in- 
dicates the 2a confidence limits. The overall best fit is obtained 
for a = O-SSIq'^ (95 per cent confidence limits). 



From the summed x data we obtain a best-fitting inner 
slope of a — 0.881q'3^, where the uncertainties are 95 per 
cent confidence limits. We conclude that our combined result 
on the inner dark matter density slope is consistent with the 
range of values predicted by CDM simulations. 



4 DISCUSSION 

4.1 Previous results on the inner dark matter 
density slopes 

Several of the clusters in the present sample have been the 
subjects of previous work that has examined the issue of 
the inner dark matter den sity slope. MS2137.3-23 53 and 
Abell 963 were studied by ISand et all (|200l l2004h using 
a combination of strong gravitational lensing data and opti- 
cal velocity dispersion measurements for the dominant clus- 
ter galaxies. These authors measured surprisingly flat inner 
slopes for MS2137.3-2353, Abell 963 (which are formally 
consistent with the findings from the Chandra data pre- 
sented in Table 3J) and four other clusters, and concluded 
that their results were inconsistent with the standard NFW 
model (a = 1.0) at 99 per cent confidence. 

Later work by Bartelmann fc Meneghettil l|2004t ) and 
iDalal fc Keetonl (|2003h cautioned that allowing for devia- 
tions from axial symmetry in the strong lensing analysis, 
using elliptical rather than spherical mass models (as mo- 
tivated by the data) will modify the constraints on the in- 
ner dark matter profil es. These authors conclu ded that the 
strong lensing data of Sand et al. 

I (|2002j . I 20041 ) remain con- 
sistent with CDM models, once such effects are taken into 
account. The X-ray results presented here are also consis- 
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tent with the standard NFW model (a = 1.0) at 95 per 
cent confidence. It is important to note that X-ray data do 
not suffer in the same way from uncertainties regarding ax- 
ial symmetry. For the regular, apparently relaxed clusters in 
the present sample, the mass models determined from the 
X-ray data under the assumption of spherical symmetry are 
unlikely to be affect ed by asphericity effects by m ore than a 
few per cent dPiffaretti et alj|2003l ; iGavazzi 2005). 

Using Chandra data for Abell 2029 (the first ~ 20ks ob- 
servation only) and parameter ized models for th e gas den- 
sity and temperature profiles, iLewis et al.l J2003) obtained 
q = 1.19±0.04 for the inner (total) mass profile. We measure 
a = 1.16^034 f° r the dark matter at 68 per cent confidence 
using 94ks of clean data. (For the total mass, a fit with a gen- 
eralized N FW model g i ves a = l.lijj^jg). Thus, the present 
study and lLewis et al. (|2003T ) obtain similar best-fit results, 
although the statistical uncertainties reported here are sig- 
nificantly larger, despite being based on more data. In part, 
this highlights the effects that priors, in the form of parame- 
terized models for the gas density and temperature profiles, 
can have on the ana lysis. IZappacosta et all (|2006l . see also 
iBuote fc LewisfeoO^ ) present a detailed study of the nearby, 
intermediate temperature cluster Abell 2589, for which they 
also model the dark matter, X-ray emitting gas and stellar 
mass associated with the cD galaxy separately. They con- 
clude that the standard NFW model with a — 1 provides a 
good description of the dark matter and total mass distri- 
but ions in the cluster. 

lArabadiis fc Bautj (2004) studied Chandra data for 
Abell 1835, Abell 2029, Abell 2204, Zwicky 3146 and MS 
2137.3-2353. Fig. 4 of that paper indicates best-fitting values 
of a ~ 1.85 for Abell 2029, a ~ 1.8 for Abell 2204, a ~ 0.9 
for Abell 1835, a ~ 1.7 for Zwicky 3146 and a = 1.6 ± 0.2 
for MS 2137.3-2353 (with 68 per cent confidence limits of 
~ 10 — 20 per cent) . In general, we find shallo wer inner dark 
matter slopes than lArabadiis fc Bautzl (|2004T i. Our statisti- 
cal uncertainties are also la rger in some cases. 

IVoigt fc Fabianl (|2006h present results on the inner 
slopes for 12 clusters, 8 of which are in common with the 
present study. These authors use a similar non-parametric 
spectral deprojection technique to measure the deprojected 
X-ray temperatures in the clusters. Their results on the cen- 
tral slope are in general agreement with those presented 
here. 



M/Lv value). We have adopted these parameters as a tem- 
plate to approximate the mass contribution from stars in 
the central galaxies of all clusters in the sample, employing 
either a Jaffe model or de Vaucouleurs model. 

In principle, we might expect the results on the inner 
dark matter density slopes to be sensitive to the choice of 
parameters used to describe the central stellar mass com- 
ponents. Here, the main systematic uncertainty lies with 
the assumption that the central stellar mass distribution 
in MS2137. 3-2353 provides a reasonable model for other 
galaxy clusters in the sample. The assumption of a con- 
stant central stellar mass is well motivated; K-band obser- 
vatio ns of dominant clus ter galaxies in X-ray luminous clus- 
ters (|Brough et al.| [2002) show little variation in total stellar 
mass over the redshift range < z < 1, with scatter at the 
level of 20 — 30 per cent. 

To estimate of the effect of departures from a con- 
stant central stellar mass on the measured inner dark matter 
slopes, we have re-analysed the Chandra data for MS2137.3- 
2353 varying the M/Lv ratio over the range M/Lv = — 6 
(and thereby changing the central stellar mass by ~ ±100 
per cent). In each case, we have re-determined the con- 
straints on the inner slope a. The results are shown in Fig. [4] 
Note that MS2137.3-2353 is one of the least massive clusters 
in the sample, and so changes in a with central stellar mass 
are likely to be larger than for most other clusters in the sam- 
ple. For M/Lv = 1.5, we measure a — I.IOIqIso a ^ 68 per 
cent confidence. For M/Lv = 4, we measure a = 0.94lQ4g. 
The tendency is for higher M/Lv ratios to give slightly 
shallower inner dark matter slopes. These results can be 
compared to the best-fitting result for MS2137.3-2353 with 
M/Lv = 2.75 of a = l.OOt^g. We conclude that the re- 
sults on the inner dark matter slope for MS2137. 3-2353 are 
robust against changes in the central stellar mass by ±50 
per cent. 

Finally, we note that extending the analysis to model 
the mass components associated with stars external to cen- 
tral galaxies separately should have a negligible effect on the 
results. In total, stars co ntribute only ~ 2 per cent of the 
total cluster mass (e.g., iLin fc Mohrll2004 iFukugita et al.l 
Il998t) 

In conclusion, the results presented here are robust 
against systematic uncertainties associated with measuring 
the stellar mass contribution in the clusters. 



4.2 On the robustness of the central galaxy model 

ISand et al l (|2002l ) present a detailed analysis of the optical 
pro perties of the dominant galaxy in MS2137. 3-2353 (see 
also Gavazzi 2005). These authors fit a de Vaucouleurs pro- 
file to the galaxy surface brightness, yielding an effective ra- 
dius R c = 5.02 ± 0.50 arcsec and a total V-band luminosity 
Ly = 4.2 x 10 11 Lq. Using the redshift-de pendent V-band 
mass-to-light ratio of (|Fukugita et al.|[l998h and assuming a 
galaxy formation redshift ~ 2, we obtain a V-band mass- 
to-light ratio for the dominant galaxy in MS2137. 3-2353 of 
M/Lv = 2.75 and a total stellar mass of 1.14 x 10 12 M Q . 
(The stellar mass dominates the tot al mass within 10 kpc of 
the center of the cluster. Note that I Sand et alj |2002) mea- 
sure an optical velocity dispersion for the dominant galaxy 
in MS2137.3-2353 of a ~ 275 km/s, f rom which they in- 
fer M/Lv ~ 3.1; see also iGavazzil 12005] for a slightly lower 



4.3 Residual systematic uncertainties 

The clusters studied here are the most massive, dynamically 
relaxed clusters known and are, therefore, the systems for 
which the assumption of hydrostatic equilibrium should be 
most robust. Although, as discussed above, geometric effects 
and uncertainties associated with separating the dark mat- 
ter and baryonic matter components are unlikely to impact 
on the results significantly, some systematic uncertainties 
remain. 

In the first case, it remains possible that small levels of 
non-thermal pressure support due to e.g. gas motions, cos- 
mic rays or magnetic fields could be present in the X-ray 
emitting gas and bias the measured masses low. However, 
the relaxed nature of the clusters argues that large, unac- 
counted for, bulk and/or turbulent motions are unlikely to 
be present and that the mass measurements should be re- 
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Figure 4. Goodness of fit for the inner slopes a for the galaxy 
cluster MS2137. 3-2353. The mass model includes a Jaffe compo- 
nent for the luminous matter of the cD galaxy and a generalized 
NFW component for the dark matter halo. 



covered to better than 10-20 per cent accuracy |Nagai et al.l 
120071 ; |R asia et all [2006). Gal axy clusters are know n to con- 
tain magnetic fields (e.g., ICarilli fc Tavlod 120021 ) The is- 
sue of magnetic pressure support and its effect on X-ray 
mass measurements has been studied bv lDolag fc Schindleij 
(2000) who showed that for relaxed clusters, magnetic pres- 
sure support is unlikely to bias mass measurements signif- 
icantly, even in the central regions. For individual clusters, 
effects no larger than 10-20 per cent are expected. The sim- 
ilarity of the mass results obtained using method 1, where 
we measure the total mass and method 2, where we model 
the dark and baryonic components separately, also argues 
that any non-thermal pressure component, at the ~ 10 — 15 
per cent level, distributed in a similar manner to the dark 
or baryonic matter, is unlikely to have a significant effect 
on the conclusions. A program to determine the maximum 
size of non-thermal pressure in the clusters, using a combi- 
nation of Chandra X-ray data and wide field weak gravita- 
tional lensing observations, is underway (Donovan et al. , in 
preparation). 

It should also be recognized that although we quote re- 
sults on the mass-concentration relation appropriate for the 
virial radii in the clusters, thereby allowing an easy compar- 
ison with the predictions from simulations, the X-ray data 
only extend to about r2soo or approximately a quarter to 
a third of the virial radius in most clusters. Some system- 
atic uncertainty is associated with extrapolating the allowed 
range of NFW mass models to these larger radii. 

Finally, we reiterate that the simulations used to pre- 
dict the mas s-concentration relation and central dark matter 
slopes (e.g., Navarro et al.lll995l . Il997l : iDiemand et al 1 l2004 



Bower et al.ll2006l ) proper- 



X-ray emitting gas and stars, and realistic baryonic physics 
(cooling, star formation, AGN feedback). For the most accu- 
rate comparison with the data presented here, such simula- 
tions should contain a sufficiently large number of massive, 
relaxed clusters and be normalized to match the observed 
X-ray (e.g., te mperature profiles, virial relations, X-ray gas 
mass fractio n; lAllen et al. l200ll |2004 IVikhlinin et aLll2006l ) 
and optical (|Croton et al. I 20061 ; 
ties. 



5 SUMMARY 

We have used the Chandra X-ray Observatory to measure 
the dark matter and total mass profiles for a sample of 34 
massive, dynamically relaxed galaxy clusters. Our analysis 
has employed a non-parametric, spherical deprojection tech- 
nique that minimizes the need for priors associated with pa- 
rameterized models for the X-ray gas density and/or temper- 
ature profiles. This allows a direct assessment of the good- 
ness of fit to the Chandra data provided by a variety of 
simple mass models as well as an accurate determination of 
statistical uncertainties on fit parameters. 

We have shown that the NFW model, which is moti- 
vated by CDM simulations, provides a good description of 
the total mass and dark matter distributions in the major- 
ity of clusters. In contrast, the singular isothermal sphere 
model can, in almost every case, be firmly ruled out. Com- 
bining the results for all clusters for which the NFW model 
provides an acceptable description of the dark matter pro- 
files, we obtain a best-fitting result on the inner slope of the 
dark matter density profile in the clusters, a = Q.88l '*i (68 
per cent confidence limits). 

We observe a well-defined mass-concentration relation 
for the clusters with an intrinsic scatter in good agreement 
with the predictions from simulations. The slope of the mass- 
concentration relation, c oc M° ir /(1 + z) b with a = —0.45 ± 
0.12 at 95 per cent confidence, is significantly steeper than 
the value of a ~ —0.1 predicted by CDM simulations for 
lower mass halos. The redshift evolution, b — 0.71 ± 0.52 at 
95 per cent confidence, is consistent with the value b ~ 1 

predicted by those simulations. 

After this work was preprinted on astro-ph. lBuote et al.l 

(2006) preprinted a paper in which they discuss the mass- 
concentration relation for 39 systems with masses in the 
range 6 x 1O 12 M to 2 x 10 15 M Q . The results for higher 
mass systems are drawn fro m the previous studi es of 
iPointecouteau et alj (|2005l ) and IVikhlinin et~aTl (|2006l ) and 
are in good overall agreement with the present work. 



Sh aw et al1l2006l ) are CDM only. In detail, these predictions 
may be modified by future numerical work that includes the 
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